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Abstract 

Gaussian QCD sum-rules are ideally suited to the study of mixed states of gluonium (glueballs) and quark (qq) 
mesons because of their capability to resolve widely-separated states of comparable strength. The Gaussian QCD 
sum-rules (GSRs) for all possible two-point correlation functions of gluonic and non-strange 1 = quark scalar 
{J^'-^ = 0+"*") currents are calculated and analyzed. For the non-diagonal sum-rule of gluonic and qq currents we 
show that perturbative contributions are chirally-suppressed compared to non-perturbative effects of the quark 
condensate and instantons, implying that the mixing of quark mesons and gluonium is of non-perturbative origin. 
The independent predictions of the masses and relative coupling strengths from the non-diagonal and the two 
diagonal GSRs are remarkably consistent with a scenario of strongly-mixed states with masses of approximately 
1 GeV and 1.4 GeV. The mixing is nearly maximal with the heavier mixed state having a slightly larger gluonium 
content than the lighter state. 

1 Introduction 

The interpretation of the nature of the lightest scalar mesons is one of the most fascinating problems in hadronic 
physics. The plethora of scalar {J^'~^ — 0++, / = 0) states below 2 GeV 1 cannot be described by a simple qq nonet, 
a situation indicative of exotic states such as gluonium (glueballs) or multi-quark {qqqq) states amongst the known 
scalar mesons. In the case of gluonium states, analyses based on chiral Lagrangians 2, 3 suggest that the /o(1500) 
and /o(1710) are mainly gluonium states with a small gluonium component of the /o(980) Similarly, QCD 
Laplace sum-rules based upon (diagonal) correlation functions of cither gluonic or qq currents find that admixtures 
of gluonium and qq exist with masses of approximately 1 GeV and 1.6 GeV [31[51[5], a conclusion that is also upheld 
by studies based on Gaussian QCD sum-rules [71 [51 [S]- Phenomcnological analyses present a scenario of mixing 
between a 1 GeV glueball and the /o(980), /o(1500) states of a qq nonet [10]. Lattice QCD calculations lead to a 
scalar gluonium state of approximately 1.6 GeV with quenched quarks [llj . However, with dynamical quarks the 
mixing with qq states appears to be very strong, driving the mass of the lightest "flavour-singlet" meso down 
toward 1 GeV with tentative identification of an excited state on the order of 1.5 GeV [T^ . 

The results of these different approaches suggest a consistent scenario where the mixing of qq and gluonium 
is manifested in the scalar hadronic spectrum as a lighter state on the order of 1 GeV and a heavier state on the 
order of 1.5 GeV. From both the QCD sum-rule and lattice perspectives, this implies that states exist which couple 
to both gluonic currents and to quark currents, and hence the non-diagonal correlation function between qq and 
gluonic currents must be non-trivial. Surprisingly, the study of such non-diagonal correlation functions has been 
largely ignored apart from the use of low-energy theorems [TBj ; to our knowledge Ref. [TJ is the only study of non- 
diagonal gluonic-gg correlation functions and its detailed analysis is limited to the pseudoscalar case. Analysis of the 
non-diagonal correlation functions is thus essential for a complete understanding of the scalar gluonium-q'g system. 

Gaussian sum-rules have been shown to be sensitive to the hadronic spectral functions over a broad energy range, 
and analysis techniques have been developed to exploit this dependence to determine how resonance strength is 
distributed in the spectral function [3 [HI [S] ■ Thus Gaussian sum-rules are well-suited to situations such as qq- 
gluonium mixing where multiple hadronic states could contribute to a correlation function. 
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The purpose of this paper is to demonstrate that all possible Gaussian sum-rules for two-point correlation functions 
amongst gluonic and qq currents {i.e. diagonal gluonic, diagonal qq, non-diagonal gluonic-qq) lead to a consistent 
scenario for strong (near-maximal) mixing of gluonium and qq scalar mesons. The formulation and analysis of 
Gaussian sum-rules is reviewed in Section [21 and previous results for the diagonal gluonic and qq Gaussian sum-rules 
are presented. In Section [3] the leading-order perturbative, QCD condensate, and instanton contributions to the 
non-diagonal correlation function of qq and gluonic currents are calculated along with the associated Gaussian QCD 
sum-rule. The analysis of the Gaussian sum-rules and the pattern of state coupling mixing is then presented in 
Section H 



2 Review of Gaussian Sum-Rules 

The simplest Gaussian QCD sum-rule (GSR) has the form [15] 

1 



Go(s,t) 
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'-{t-sf 


/exp 
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to 





TT 



p{t)dt , r>0 (1) 



and relates the QCD prediction Gq (s, r) to an integral of its associated hadronic spectral function p{t). The smearing 
of the spectral function by the Gaussian kernel peaked at t = s through the (approximate) region s — 2^ < t < 
s -|- 2-y/r provides a clear conceptual implementation of quark-hadron duality. The width of this duality interval is 
constrained by QCD since renormalization-group improvement of the QCD (left-hand) side of ([1]) results in identifying 
the renormalization scale v through — ^/r [3 [15] ; therefore it is not possible to achieve the formal r — > limit 
where complete knowledge of the spectral function could be obtained through 

lim Go (s, r) = -p {§) , s > to . (2) 

The variable s in ([T]), on the other hand, is unconstrained by QCD, and so the s dependence of Gq (s, r) can be used 
to probe the behaviour of the smeared spectral function, and hence the essential features of p(t) can be extracted. 

An interesting feature of the GSR ([T|) is its ability to study excited and ground states with similar sensitivity. For 
example, as s passes through t values corresponding to resonance peaks, the Gaussian kernel reaches its maximum 
value. Thus any features of the spectral function strong enough to be isolated from the continuum will be revealed 
through the GSR. In this regard, GSRs should be contrasted with Laplace sum-rules 

oo 

i?(A^) = i y'exp(^-^)p(t)dt , (3) 

to 

which exponentially suppress excited states in comparison to the ground state. 

Gaussian sum-rules are based on QCD correlation functions 11 [Q^) of renormalized composite operators Ji{x) 
and J2(a;jl 

n (Q2) ^ijd^x e'«- (0|r [Ji(x) J2(0)] |0) , ^ ~q' . (4) 

By allowing the possibility that Ji ^ J2, our formulation encompasses both diagonal (Ji — J2) and non-diagonal 
correlation functions. The correlation function ((4|) will satisfy dispersion relations appropriate to the asymptotic 
form of the correlator in question. For example, the diagonal scalar gluonic correlation function 11 (Q^) satisfies the 
following dispersion relation with three subtraction constants 

C30 

n (q2) - n(o) - Q2n'(o) - \q'i^'{o) = -^J -^t^qt^ dt . (5) 

to 

In general, the quantity p{t) is the spectral function appropriate to the quantum numbers of the given currents, and 
it should be noted that in certain situations (such as the scalar gluonic correlation function) the subtraction constant 



^The configuration-space correlation function in (|3J as used in lattice QCD will contain an exponential suppression of excited states 
similar to that occurring for Laplace sum-rules. 
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n(0) is determined by a low-energy theorem [T3]. Undetermined dispersion-relation constants and field-theoretical 
divergences are eliminated upon constructing the GSR [7^1 



Gk(s,T)^ f + - - - iA)'-n(-g + iA) 

where A: = — 1,0,1,... and where the Borel transform B is defined by 

A2/W=4t 

The dispersion relation ([5]) in conjunction with definition ^ together yield the following family of GSRs 
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(9) 



has been used to simplify the phenomenological (right-hand) side of ([5]) as well as the term on the QCD side 
proportional to 11(0). Clearly the k — —1 sum- rule can only be defined in cases where there exists an appropriate 
low-energy theorem because 11(0) is present in ([5]) when k = —1. Calculation of Gk{s,T) is achieved through an 
identity relating ([7]) to the inverse Laplace transform ^15j 



6[/(A2)]^-£-i[/(A^)] 
where, in our notation, 



(10) 



(11) 



with a chosen such that all singularities of / lie to the left of a in the complex A^-plane (see [7] for further details). 
Next, we impose a fairly general resonance(s) plus continuum model 



p{t) = phad(i) + 61 (i - So) ImnQCD(t) 



(12) 



where sq represents the onset of the QCD continuum. The continuum contribution of (jl2p to the right-hand side 
of m is 



Gr'{s,T,so) 
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' Attt 



t'' exp 



At 



and is combined with Gk (s, r) to obtain the total QCD contribution 

G?^° (s, r, so) ^ Gk {s, t) ~ G^' (S, r, so) , 
resulting in the final relation between the QCD and hadronic sides of the GSRs 



G^^° (s,r. So) -1-4- 
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/Attt 
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Original studies involving the GSRs exploited the diffusion equation 
a2Go(s,T) _ 5Go(s,t) 



dp 



dT 



(13) 



(14) 



(15) 



(16) 



^This definition is a natural generalization of that given in I15| . To recover the original GSR, we simply let fc = in ||6j. 
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which follows from ([5]) for fc = 0. In particular, when p^^^(t) [see (IT^ ] is evolved through the diffusion equation (|16p. 
it only reproduces the QCD prediction at large energies (r large) if the resonance and continuum contributions are 
balanced through the A: = member of the finite-energy sum-rule family [15j 

So 

Fk{so)^lJt'p''^''{t)dt . (17) 

to 

An additional connection between the finite-energy sum-rules (IT7|) and the GSRs can be found by integrating 
both sides of with respect to s to obtain 

oo oo 

y G^^°(s,r,so)ds + 4-in(0) = iytV"^(i)dt , (18) 

— oo to 

indicating that the finite-energy sum-rules (|17p are related to the normalization of the GSRs. Thus the information 
independent of the finite-energy sum-rule constraint following from the diffusion equation analysis |T3] is contained 
in the normalized Gaussian sum-rules (NGSRs) [8] 

^gco,^ , Gg°''('.^.».)+-i^-i^<''p(Tr)n(0) 

fQCD/ N _ / ^QCD 



Mfe^^"(T,5o)= J G^^"(s,r,so)ds , n = 0,l,2,... , (20) 

— OO 

which are related to the hadronic spectral function via 



For diagonal correlation functions the spectral function obeys a positivity constraint so the NGSR must exist. 
For non-diagonal correlators the possibility of state mixing implies that p^^'^ (t) could change sign, so it is possible 
that either M^^°(t, So) or that the denominator on the right-hand side of (PT|) could be zero. In such situations, 
the GSR HID would have to be analyzed instead of the NGSR (EH). 

We now consider the currents that will be used to probe the gluonic and qq aspects of the hadronic states. 
Refs. [H [TB] argue quite eloquently that the mixing of qq mesons and gluonium is unavoidable because of the trace 
anomaly for the energy-momentum tensor T^i, |17| 

r; = ^(]{a)Gl,G'^^-^ + [1 + 7(a)] J2 "^*^/^/ (22) 

/ 

where 

11 1 ^51 19 

/5o = ^-gn, , = (24) 

„ , , , dm 

-2m7 a = v^— . (25) 

In fact, (|22p contains two multiplicatively-renormalizable (renormalization-group invariant) composite operators miptp 
and f3G^ + 4:'ym'tjj'tp , so from a field theoretical perspective these two operators would be suitable choices for currents. 
However, with this choice the gluonic or qq nature of states coupled to the current f3G'^ + 4:'jmipip is obscured. Instead, 
we follow 14J and choose the renormalized currents 

Jg{x) ~ aG\ = a ^1 + — — ^ G\ — Aa \muu{x)u{x) + mdd{x)d{x)\ ^ + . . . (26) 

Jq{x) = rriq [u{x)u{x) + d{x)d{x)] , niq — (m„ + nid) /2. (27) 
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where the notation R refers to the renormahzed composite operator and the B denotes an expansion in terms of 
bare quantities and our convention for dimensional regularization uses D = 4 + 2e spacetime dimensions. Of course 
the form of the renormahzed operator (|26p necessarily underlies the renormalization-group invariance of the trace 
anomaly ([221) (see e.g. [18j). However, the advantage of the current Jg is that its tree-level expansion is purely 
gluonic allowing a qualitative separation of gluonic and qq degrees of freedom. The current Jq has isospin / = and 
is renormalization-group invariant. 

The diagonal correlation function for the gluonic currents Jg is defined by 



(Q') = i J d^xe'^-- {0\T[Jg{x)Jgm |0) , 



(28) 



while the diagonal correlation function of the / = (non-strange) quark currents Jq is 

n,, (Q") - i / d'^: e'«- (0|r [Jg{x)Jg{Q)] |0) , = -q\ (29) 

Although both the gluonic and (/ — 0) quark diagonal correlation functions are probes of scalar hadronic states, those 
states which have a more significant overlap with the gluonic current should predominate in (|28p. while those states 
which are dominantly of a quark {qq) nature should be more significant in (|29p . A mixed state with substantial 
gluonic and quark components {i.e. a state that overlaps with both the gluonic and quark currents) should self- 
consistently appear in an analysis of both correlation functions. In particular, prediction of mass-degenerate states 
from the QCD sum-rule analysis associated with these two correlation functions is indicative of a mixed state. 
In the scalar gluonic channel, the low-energy theorem (LET) [T3j 



ns9(0)^ Jim^n3g(Q2) 



87r 



iaG' 



(30) 



allows construction of the fc = — 1 GSR. The significance of instanton contributions in the overall consistency of 
the LET-sensitive k = —1 sum-rule and the LET-insensitive fc > sum-rules was first demonstrated for Laplace 
sum- rules [191 120] . A similar consistency is observed for the Gaussian sum-rules, but theoretical uncertainties are 
better controlled in the fc > GSR 7 , and hence this paper will focus on the fc = GSRs for the diagonal gluonic 
and quark channels. 

The QCD correlation functions (p8)) and (p9| contain perturbative, condensate, and instanton contributions. For 
the diagonal scalar gluonic case, the QCD prediction for the fc = GSR to leading order in the quark mass is 
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(31) 



The perturbative coefficients in (|3ip are given by 



ao 



659 a 
"36"^ 



247.480 



ai 



65.781- 

4 TT 



0-2 



-10.1250 f- 

VTT 



(32) 



as obtained from the three- loop MS calculation of the correlation function in the chiral limit of n/ =3 flavours [2T] . 
The condensate contributions in (j3ip involve next-to-leading order 22J contributions from the dimension four gluon 



*To leading chiral order the GSR is identical for the current II26I I or the trace-anomaly version /3G^ -I- A-yrnqq. 
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condensate (aG^) and leading order '23] contributions from gluonic condensates of dimension six and eight 



1 A ^ 
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Co = 8tt^ (- 

\ TT 
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do 



TT 



(33) 
(34) 

(35) 



The remaining term in the GSR pip represents instanton contributions obtained from single instanton and anti- 
instanton [21] {i.e., that multi-instanton effects are negligible [25j) contributions to the scalar gluonic correlator [19, 
I2Q[ I23[ I26j within the liquid instanton model [19] parameterized by the instanton size p and the instanton density 
ric. The quantities J2 and Y2 are Bessel functions in the notation of p7] . 

As a result of renormalization-group scaling of the GSRs [HI [H] , the coupling in the perturbative and condensate 
coefficients and (1551) ] is implicitly the nmning coupling at the scale v'^ ~ \Jt for n/ = 3 in the MS scheme 

a(}?) _ 1 /3ilogL 1 



/3oi (/?oi)' (/3oi)' 



log 
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/3o 



[/32 (log^L-logL-l) +;32] 

3863 



(36) 



/3i 



/32 



384 



with A jj-g « 300 MeV for three active flavours, consistent with current determinations of a(MT-) [T]. 

Analogous to the scalar gluonic case, the QCD prediction of A; = diagonal GSR for the 7 = scalar quark 
currents to leading order in the quark mass is [8] 
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(37) 



Again, renormalization-group improvement implies that both nig and a are implicitly running quantities at the scale 
= y/r as given by ([36]) and the (two-loop, n/ = 3, MS) expression 



290 1 256 log L 
729 L ~ 729 L 



log 



(38) 



where rhq is the renormalization-group invariant quark mass parameter. The perturbative contributions in (I37p are 
the 71/ = 3 two-loop results obtained from [|28 ■ , and the instanton expressions are obtained from [29] . The condensate 
contributions are leading-order results obtained from [29 , , and are defined by the quantities 



{GlOl 



gOT> + (ttsG^) , (qq) = ^{uu + dd) 



and 



u,d,s I 



(39) 



(40) 



The vacuum saturation hypothesis [29] in the SU (2) limit {uu) = (dd) = (qq) provides a reference value for (Og) 

(Gg^O^) = -Us^as {{qqf) - -/..1.8 X 10-*GeV« , (41) 



where the quantity fyg parameterizes deviations from vacuum saturation where fys — 1. 

The GSRs ([3T|) and ([37)1 exhibit some interesting qualitative features. For example, the condensate contributions 
decay exponentially with the Gaussian peak-position s, emphasizing that these contributions have a low-energy 
origin. 
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3 Non-Diagonal Correlation Function and GSR of qq and Gluonic Cur- 
rents 



The non-diagonal correlation function for quark and gluonic currents is defined by 

Hg, (Q') = i J d^xe">-- {0\T[Jg{x)Jgm |0) , = . (42) 

As with the diagonal correlators, the non-diagonal correlation function contains perturbative, QCD condensate, and 
instanton contributions. 

The leading order perturbative diagrams that contribute to Ugq (Q^) are given in Figure [TJ The first diagram is 
a two-loop calculation and corresponds to the (bare) gluonic term in (p6)) while the second is a one-loop calculation 
corresponding to the (bare) quark term arising from composite-operator renormalization. However, one can see that 
both diagrams lead to the same order in Both diagrams also have the same O (wg) leading chiral behaviour; the 
first diagram results in an O (m^) chiral suppression factor from the quark loop. 



Figure 1: Leading order perturbative diagrams for the non-diagonal correlation function. The symbol (8" denotes the 
bare current aG^ within Jg and © denotes the bare current Jq. In the second diagram, the solid square represents 
the qq term arising in the renormalization of Jg. The Feynman diagrams were drawn with JaxoDraw [30| . 



(43) 



The two- loop diagram in Figure [T] can be expressed in terms of the fundamental two-loop integrals [31 
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(46) 
(47) 



where we have converted the dimensional regularization conventions in fST] to D = A + 2e. The integral Iq is strictly 
not necessary since it is finite, but we provide it for completeness. At leading chiral order, the contribution of the 
two-loop perturbative diagram to Ugq is 



USna^m^ [-8I2 

2 2 ^ ^ 

a m 7; 



2D)p'h + {^D - 8)p2/4 + (4 - 2D)h + (2 - D)p^h] 



2 ^ 

+ a m 5— 



L + 27- 



35 



log independent terms , L = log (Q'^/4t 



(48) 



^Strange (and heavier) quarks originating from composite operator renormalization in II26I I will be suppressed by a factor of a^. 
®Note that the integral does not appear in the the classification of Ref. I31| : it has been included here for completeness. 
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The non-logarithmic terms in (j48p can be ignored since they correspond to dispersion relation substraction constants 
which are eliminated when forming the GSR. 

The L/e term in (|48p is problematic since it cannot be rcnormalized away or absorbed into a dispersion-relation 
subtraction constant. However, the leading chiral order contribution of the one-loop diagram of Figure [1] to 11^^ is 

48q2to2q2 

-I —Jl , (49) 

where we have introduced the following notation for one-loop integrals: 
f d^k 1 

^1-/7:^1:57:-^' (50) 
= I .o_^D .2.„ .^4 = / ^2n)D k\q - kf ' ^^^^ 
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(27r)^ 


fc2(q- 


A:)2 ' 
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/c)4 


d^k 
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/fc)4 ■ 



~ i (2^)^ k^{q - kY ■ ^^^^ 
Thus the contribution of the one-loop "renormalization" diagram to Ilgq is 
2 2 iQ^L 1 



-a m 



log independent terms , (53) 



which cancels the problematic L/e term from the two-loop diagram. Thus the leading-order perturbative contribu- 
tions to the non-diagonal correlation function are 

(Q') = ^^Q^ [^oi + A^L^] (54) 

3/a\2/ 35 \ . 3 
^0 = - - 

TT V TT 



This result disagrees with that presented in Ref. [2], although the overall chiral and logarithmic dependence is 
qualitatively similar. We have checked our calculational methodology by verifying the pseudoscalar results in 
and believe that ([55)1 is correct 

The quark condensate contributions to the non-diagonal correlator are easily calculated to leading chiral order 
using any of the equivalent methods for evaluating operator-product expansion coefficients [32j . As in the perturbative 
case, in principle there are the two classes of diagrams shown in Figure [2l However, the result of the first set of 
diagrams of Fig. [2] is 

mi\a^mq{qq) [{~2D - 4) Ji -f {2D - A)q^ J2 + (2 - D)q'J:i\ = --a^mq{qq) log {%\ (56) 

where logarithm-independent terms have been ignored as they correspond to dispersion relation substraction con- 
stants which are eliminated when forming the GSR. The "renormalization-induced" diagram of Fig. [2] is chirally- 
suppressed relative to (|56| so it represents subleading effects. 

It is more difficult to calculate the gluon condensate contributions, which may explain why the non-diagonal 
correlators have not previously been the subject of detailed study. The diagrams that could lead to the gluon 
condensate contributions are shown in Figure [31 However, the "renormalization-induced" diagrams are higher-order 
in a and hence are subleading. Because of infrared divergences, it is necessary to retain the quark mass until the last 
steps of the calculation and then extract the leading chiral behaviour. Using plane-wave methods with rriu = rrid = m 
we find 

327raTO^ .-, r n 

2i—^{aG') [m^ (8 - 4D)Ki - 2DK2 + {D^ - 4D + 8)^3 , (57) 



''Ref. |14j only provides the final result for the scalar case; Eq. I|48ll presents the results in a form similar to the pseudoscalar results 
contained in 1141 . 
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Figure 2: Leading order quark condensate diagrams for the non-diagonal correlation function. The solid circles on 
the quark lines denote insertion of plane-wave states or coordinate space vacuum expectation values for evaluation 
of the operator-product expansion coefficients. All other notations are identical to Fig. [T] 



where 



1 



r(2 + e)r(-e) 
r(2 + e) ' 



log 



^1 -4m2/q2 _ 1 



(27r)^ (^^ - "i^) [(? - - m2 
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^1 - 4m2/g2 _ 1 



(58) 
(59) 

(60) 



The divergences in and K2 cancel, leaving a finite result as required given the m? jq^ pre- factor. The logarithmic 
correction from Ki is seen to be subleading compared with the logarithmic correction from /^a. Thus the leading- 
chiral gluon condensate contribution to the non-diagonal correlator is 



2am^ , „, 
^ aG^ 



7rQ2 



3 - log 



(61) 







Figure 3: Leading order gluon condensate diagrams for the non-diagonal correlation function. The solid circles on 
the gluon lines denote insertion of plane-wave states for evaluation of the operator-product expansion coefficients. 
All other notations are identical to Fig. [1] 



An additional subtlety must be addressed before we can accept ([61]) as definitive. It is known that the naively- 
calculated (aG^) results may differ from the actual (aG^) contributions when the quark mass is retained to obtain 
subleading mass corrections in the gluon condensate contributions j33j . The resolution of this problem is related to 
operator mixing; the coefficient of the operator G^ receives a correction of 



12to 



-Cc 



(62) 
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where Cqq is the coefficient of the quark condensate {qq) [53]. In our case, Cqq ~ ma^ which provides a ~ 
correction to the coefficient of , representing a subleading correction compared with (|6ip . Hence we conclude that 
no further calculation is needed to obtain the leading gluon condensate corrections. 

Finally, the result of our calculation for the single instanton [53] contributions to the non-diagonal correlator in 
the dilute instanton hquid model |19j are: 



(63) 



where Kn is a modified Bessel function in the conventions of [27] . 

Combining Eqs. (f55|) . ([56]) . (|6T|) . and ((63|) we have the leading-order chiral and a contributions to the non-diagonal 
correlation function of gluonic and (non-strange) / = quark currents: 



Hg, (Q') =m2g2 [AoL + AiL^] +m,(g<7)Colog 



Bo + Bi log 



Co = -8^ - 

V TT 

Bo = 6- , Bi 

TT 



(64) 

(65) 
(66) 



along with Aq, Ai given in ([55]) . From the correlation function, the k = GSR can be calculated as outlined in 
Section H 



I Attt 



dt exp 



At 



t[Ao + 2Ai log 



Coniqiqq) 



Bim^giaG 




1 



log 



Attt 



cxp 



to; 



s 



exp 



'4r j 



(67) 



tVt Ji (p^ft) Y2 (p^ft) + J2 (pVt) Yi (pVt 



where in practice, the limit r/ is implemented numerically with values of r/ < 10 ^ in GeV units. The normalized 
GSR is defined by 



nI,'"^ (s,r,so) = 



Mlf'^ (r,.o) 



M(^^)(r,so)= / Gr{s,T,so)d 



{a<i)i 



(68) 



In Eqs. (j67p TOg and a are implicitly the leading-order versions of the running quantities (j36|) and (j38p evaluated at 
the scale = ^/t. 

Since a low-energy theorem exists for the non-diagonal correlator |13j 



^39(0) —mq{qq) 



(69) 



the LET-sensitive k — ~1 GSR is also relevant. Again using the methods outlined in Section [2] the corresponding 
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results for the k = 



-1 GSR are 



/4ttt 



dt exp 



Comq{qq) lim 



2v3nc7rmqap 



'4ttt 




exp 



exp 




(70) 



Vt Ji (p^ft] Y2 (p^ft] + J2 ipVt] Yi (pVt 



Note that the last term in (j70p has a functional dependence identical to the LET term in (jT9|) . and hence there is 
an LET-like instanton contribution for the non-diagonal GSR similar to the diagonal gluonic case [7j. The k = —1 
NGSR is defined by 



G^ff (j,r,.o) + ^exp(^)n,,(0) 



cLf (s,r, so)ds. 



(71) 



The non-diagonal correlator and its associated GSR have quite distinct chiral behaviour compared with the di- 
agonal correlators. In the diagonal case, the perturbative, condensate, and instanton corrections all appear with 
identical powers of the quark mass. However, in the off-diagonal case the perturbative and giuon condensate correc- 
tions are chirally-suppressed compared with the quark condensate and instanton terms. The physical implications 
of this result will be discussed below. 



4 Analysis of Gaussian Sum-Rules for qq and Gluonic Currents 

The general strategy for analysis of NGSRs is fitting the QCD prediction with a parameterized model for p^^'^{t) 
in (I2ip . Correlation functions of vector and axial- vector qq currents can be directly related to experimental data 
(e.g., R{s)), but in the case of gluonium there is no direct connection with experimental observables. The narrow 
resonance approximation is the most common choice made for Laplace sum-rule analyses of gluonium, with either a 
single (narrow) resonance [51 [531 [311 [31] to examine the dominant gluonic state or two (narrow) resonances [H [SI [33] 
to explore the possibility of gg-gluonium mixtures. Laplace sum-rule gluonium analyses which go beyond the narrow 
width approximation include a single Breit-Wigner resonance skewed by kinematic factors |37| , and an interpolation 
between the LET and continuum behaviour [16] . Finite-energy sum-rule analyses of scalar gluonium that incorporate 
resonance widths include step functions [38] and Breit-Wigner resonances [39^ with kinematic skewing^ 

GSR analyses of gluonium have employed single and double narrow resonance models in addition to a variety of 
models that incorporate resonance widths [71[S|. However, since the gaussian kernel in ([T]) has a QCD-limited width 
of > 2 GeV^ , one would expect that there is insufficient resolution to observe resonance width effects and hence 
the narrow resonance may be an appropriate approximation. Such an approximation is clearly desirable for fitting 
the QGD prediction to the resonance model in (|2ip for multiple resonances. 

To explore the effect of the narrow width approximation, we consider the following models for a resonance of 
mass to: a narrow resonance 

V =/''5(t-TO2) , (72) 

a structureless square pulse of width 2TOr 
1 

-p^P(t) = -1^ \e it - + TOP) -eit-m''- TOP)] , (73) 

TT 2to1 

*Ref. 1391 also uses the Gaussian sum-rule diffusion equation (II6I I analysis to constrain the QCD continuum. As discussed in Section[2l 
our approach based on NGSRs provides information that is independent of this duality constraint. 
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and a gaussian shape 



272 



(74) 



Tbwi for the gaussian 

-\/2 log(2)7/m. The contribution of these models to the quantity appearing on the 



For the square pulse the quantity F has a direct analogy with the Breit-Wigner width F 
shape the relation is Tbw 
right-hand side of (PT|) 



1 rOC ,h. 

7= L rexp 

47rr J to ^ 



At 



p'^^d(t) dt 



can be calculated in closed form [7]: 
1 



(75) 




exp 



4t 



, s — + 77iF \ r f s — ra^ — mF 
erl — = — erl 



2^ 



1 + erf 



27r^72 + 2t 



exp 



2^/? 

(j — TO^)^ 

"2(72 + 2r) 



(76) 
(77) 

(78) 



As argued above, we see that in the limits 2-y/r ^ mF and 2\/r ^ 7 the models ((77)) and (|78p incorporating a 
resonance width will be indistinguishable from the narrow-resonance case. Fig. 2] compares (rf5)) -([75 | along with a 
numerical approximation for a Breit-Wigner shape when m = 1 GeV and Tbw = 300 MeV. The Figure illustrates that 
the resolution of the gaussian kernel is the dominant effect; resonance width effects correspond to a slight broadening 
compared to the narrow resonance case. Furthermore, resonance width effects would not alter the fitted value of 
the resonance mass. This has been explored by finding the mass M in the various models with Tbw = 300 MeV 
that provides the least-squares fit to a m = 1 GeV narrow resonance. In all cases (square pulse, gaussian shape, 
Breit-Wigner) the fitted mass M coincides with 1 GeV. 
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Figure 4: Comparison between Nq(s,t) for m — 1 GeV in the narrow resonance model and the square pulse, 
gaussian, and Breit-Wigner models with Tbw = 300 MeV. The upper set of curves are for t = 2 GeV* and the 
bottom set of curves is for r = 4GeV^. 

Our observation that the mass would be unaltered when fitting the various models to a QCD prediction can be 
understood by noting that the peaks are aligned in Fig. |3) However, if one introduces a kinematic distortion of 
into Eqs. ((73)) - ([74l) as in Refs. [13l [37l [38l [39j then the peak position will be shifted above t — rn^. For example, the 
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kinematically-skewed generalizations of Eqs. (|73p - ((7i)) are 



P^^P(i) ^ /sip 



TT 

1 
TT 



2m3r (m2 + ^ 



[e {t ~ + rat) - 6 {t - - mV)] , 



3 

2\2 



(79) 
(80) 



which lead to 

7Vo"P(s, r) =2t (s + to^ _ ^r) exp 



4^ 



ENq^{s,t) = 4 exp 



+ 2\/7r exp 



[s^ + 2r] l^erf 
5^7^ + 2Tm'' 



s — m + mF 



erf 



— 2t (s + + mF) exp 
s — m? — mF 



(s — TO — toF) 
4^ 



2^ 



(81) 



4x7^ 

(g-TO^)^ 

'2(72 + 2t) 
1 + erf 



(1^7"' + 4.s72mV + 4m''T2 + 2T"f^ + 4x^7^) 

(57^ + 2rTO2 



(82) 



where 
E 



27r(TO'* + Y) 



1 + erf 



/ to2 \ 



+ 2to 7 exp 



TO 



7^(^2 + 2r)5/2 



(83) 



In fact, the skewed square pulse is the model used in Ref. [38] for a finite-energy sum- rule analysis. If we now fit 
Eqs. (|8ip and ([5^ with Tbw = 300 MeV to find the mass M that provides the least-squares fit to a 1 GeV narrow 
resonance, M is found to be slightly less than 1 GeV (0.93, GeV and 0.98 GeV respectively). Fig. [5] shows the best 
fits; the narrow resonance still provides an excellent approximation in the presence of a kinematic skewing factor 
once the mass shift is taken into account. 
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Figure 5: Gomparison between iVo (s, r) for to = 1 GeV in the skewed square pulse and gaussian models with 
Fbvk = 300 MeV along with the best fit of a single narrow resonance to each model. The upper set of curves are for 
T — 2 GeV* and the bottom set of curves is for r = 4 GeV^. 

Thus we conclude that for purposes of a least-squares fitting methodology, the narrow-width model provides a 
good approximation to a variety of resonance models that incorporate resonance widths. However, there exists the 
possibility that the fitted mass resulting from fitting the narrow-width model may over-estimate the actual resonance 
mass by no more than 10%, i.e., the narrow- width approximation provides an upper bound on the resonance mass in 
more complicated models [40j. 
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In summary, the methodology we will employ to study qg-gluonium states including the possibility of mixing, is 
to fit the QCD prediction to the double narrow resonance model for the NGSRs in the two diagonal {qq, gg) cases 
and the non-diagonal case. In general, this methodology will take the form 



iphad(^) = fl5{t - ml) + mt - ml) , 

TT 



1 



ri exp 



(s — 



r2 cxp 



[s - ml) 



2\2-\ 



4t 



ri 



ft + fi 



r2 



fi 



f2 , ,2 ' '"1 + '"2 = 1 



ft + fi 



(84) 
(85) 
(86) 



where /i, /2 denote the couplings of the resonances to the currents under consideration. A least squares fit of the 
s, T dependence on the two sides of (|85p therefore involves minimization with respect to four parameters: the QCD 
continuum sq, the two resonance masses (mi, m2), and either ri or r2 because of the ri + r2 = 1 normalization 
constraint. The three NGSRs [see Eqs. ((3T|) . ([57]) . ([57)1 ] involving all possible combinations of the qq and gluonic 
currents will be analyzed separately. This independent analysis of each channel is essential because a signal of states 
that are t^g-gluonic mixtures would be the self-consistent appearance of states of the same mass in all three cases 
since such mixtures would necessarily couple to both the gluonic and qq currents. 

The QCD input parameters appearing within the NGSRs will now be specified. For the giuon condensate we 
employ the (central) value from [ITj 



{aC^^G"'"') ^ {aG^) (0.07 ± 0.01) GeV^ , 
and the quark condensate is determined by the PCAC relation 

mg{qq) = , U = 93MeV. 

The dimension six gluon condensate ([55)1 can be related to the value of {aG'^) using instanton techniques (see 
(Oe) = (0.27GeV^)(aG^) . 



(87) 



(88) 



(89) 



Further, by invoking vacuum saturation in conjunction with the heavy quark expansion, Ref. |42j has also related 
the dimension eight gluon condensate ([M)) to (aG^) through 



iOs)^-{{aG')y 



(90) 



We use f vs = 1.5 as a central value to accommodate the observed deviations from vacuum saturation in the dimension 
six quark condensates (|4ip ^43|. In addition, the dilute instanton liquid (DIL) model [19] parameters (which have an 
estimated uncertainty of about 15%) 



8.0 X lO""* GeV^ 



(16 



(91) 



will be employed. The NGSRs for the diagonal correlators do not require knowledge of the quark masses; in the 
diagonal gluonic case this occurs because the leading chiral behaviour is independent of the quark masses, while 
the diagonal qq case is proportional to and hence the quark mass dependence cancels when forming the NGSR. 
However, the non-diagonal case has terms of differing chiral order and therefore requires input of the quark mass. 
Unfortunately, is not known very accurately; we will use the Particle Data Group range for the 2 GeV MS mass 

m- 

2.5 MeV < mq{2 GcV) < 5.5 MeV . (92) 

The implications of the large uncertainty in within our analysis will be discussed in more detail below. 

We begin with the diagonal gluonic-gluonic case. As argued in Ref. [7] the LET-sensitive fc = — 1 GSR shows 
a strong dependence of the QCD input parameters (particularly the instanton parameters and gluon condensate) 
compared to the fc = GSR. Our current analysis will therefore focus on the fc = case@ In Ref. [7] numerically- 
efficient techniques involving the t dependence of the s peak position were developed to extract the resonance 



^The qualitative consistency between the k 
efEects are an essential feature of these works. 



and fc = — 1 cases has been demonstrated for Laplace sum-rules |19l I20| : instanton 
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and QCD continuum parameters. The numerical efficiency gained from these techniques allowed exploration of 
a variety of resonance models, including the double narrow resonance model which resulted in the central values 
So = 2.3GeV2, mi = 0.98 GeV, = 1.4 GeV, r'f<'^ = 0.28, and ^^^^ = 0.72 [HIS]. As an independent check of 
these diagonal gluonic results for this paper, we have performed a least-squares fit between the two sides of (|55|) 
in the range -4GeV2 < s < 8GeV^ and 2 GeV* < t < 4 GeV resulting in sq = 2.30 GeV^, toi = 0.951 GeV, 
m2 = 1-41 GeV, r'f^^ — 0.303, and r\^^^ = 0.697. The close agreement between the predicted resonance parameters 
and QGD continuum in these two different approaches validates the methodologies. As illustrated in Fig. [51 the 
resulting fit between the QCD prediction and double narrow-resonance model is excellent; there is no suggestion of 
resonance-width effects that would lead to differences between the QCD prediction and the phenomenological model 
on a scale similar to that of Fig. |4l The uncertainties associated with the QCD input parameters have been found 
to be 10% for r^^^^ and 0.2 GeV for the masses, with a correlated effect in the masses that lead to a relatively stable 
mass splitting mi — m2 ~ 0.5 GeV with an uncertainty of 0.03 GeV [7]. 
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Figure 6: Comparison of the theoretical prediction for the normalized GSR -/Vq^^-* (s, r, sq) with the double narrow 
resonance phenomenological model for the optimized value of the continuum sq and resonance parameters (mi , m2 , 
ri ) . The upper set of curves are for r = 2 GeV^ and the bottom set of curves is for r = 4 GeV^. The phenomenological 
and QCD predictions cannot be distinguished in the plots because they overlap almost completely. 

The diagonal qq-qq [I = 0) case does not have the complexity of a LET-sensitive sum-rule, so fc = is the 
lowest-weighted GSR that can be analyzed. The central values obtained from the s peak-position analysis are 
So = 2.6GeV2, mi = 0.97 GeV, mj = 1.4 GeV, rj"^"^' = 0.63, and 4''^ = 0.37 [HIS]- As in the diagonal gluonic 
case, we have performed a least-squares fit between the two sides of (|85p in the range —4 GeV^ < s < 8 GeV^ 
and 2GeV4 < r < 4 GeV resuhing in so = 2.348 GeV^, mi = 0.943 GeV, ms = 1.41 GeV, r(''^ ^ 0.593, and 
r^'''^^ = 0.407. Fig. [7] again shows excellent agreement of the resulting fit between the QCD prediction and double 
narrow-resonance model. Once again, there is no suggestion of anomalous differences between the QCD prediction and 
the phenomenological model indicative of substantial resonance-width effects. The uncertainty associated with r^'^^ 
is 10%, comparable to the diagonal gluonic case, but the masses show less dependence on the QCD input parameters 
than the diagonal gluonic case with uncertainties of 0.05 GeV. The uncertainty in the masses is correlated as in the 
gluonic case; the mass difference of approximately 0.5 GeV between the two states is relatively stable. 

The remarkable agreement between the resonance masses resulting from the diagonal gluonic and diagonal qq 
NGSRs suggests the existence of states with masses of « 1.0 GeV and « 1.4 GeV which are gluonium-gg mixtures, 
with the heavier state being slightly more gluonic because of its stronger coupling in the gluonic channel and a 
weaker coupling in the qq channel. The results also indicate that the mixing is rather strong because ti and r2 are 
not appreciably different, and hence the non-diagonal correlator must also contain clear signals of this strong mixing 
to validate this scenario. 

Before proceeding with a detailed analysis of the non-diagonal correlator, we consider the approximate scales 
associated with the couplings of the resonances to the gluonic and scalar currents. The perturbative corrections in 
the diagonal correlators (PT]) and (I37p imply that 

fl ^ Jl ^ mlM\ (93) 
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Figure 7: Comparison of the theoretical prediction for the normahzed GSR N^f^^ (s,t, sq) with the double narrow 
resonance phenomenological model for the optimized value of the continuum sq and resonance parameters (mi , 7712 , 
ri ) . The upper set of curves are for r = 2 GeV^ and the bottom set of curves is for r = 4 GeV^. The phenomenological 
and QCD predictions cannot be distinguished in the plots because they overlap almost completely. 



where fg and fq respectively denote the resonance couplings to the gluonic and qq currents and M is a characteristic 
sum-rule scale of order M ~ 1 GeV. In the simplest single-angle mixing scenario, the non-diagonal correlator 
will be proportional to /g/^ sin 20 where 9 is the mixing angle. The perturbative corrections to the non-diagonal 
correlator (|67p then imply 



/J,sin20^m^(^)' 
Combining (l93l) and ([94| then leads to a chirally-suppressed mixing angle 
sin 26* < 1 . 

TT M 



(94) 



(95) 



The apparent contradiction between the strong mixing found in the GSRs for the diagonal correlators and the 
basic perturbative scales in the non-diagonal GSR is resolved by a detailed analysis of the non-diagonal case. We 
first define the leading O (mg) chiral terms in (l67|) as 



Xo^''^ (s,T,so) = -Comqiqq) 



/SncTTITlqap 




dt exp 



it- 



s) 



At 



dt exp 



4t 



(96) 



tVt [ji (pVt) Y2 i^pVPj + J2 i^pVt) Yi (^pVt 



Fig.[5]demonstrates that these leading chiral terms are actually the dominant contribution to the non-diagonal NGSR, 
avoiding the chiral suppression occurring in (j95p . We thus have the intriguing result that the underlying mixing mech- 
anism is fundamentally non-perturbative, i.e., perturbative analyses do not provide the essential phenomenological 
scales. 

Because of the PC AC relation (f88| . the quark condensate term in ([96|) is independent of iriq, so in principle iVg^^^ 

could be strongly dependent on ruq. Fig. [9] shows that this is not the case; the chiral contributions to N^^''"' are 
relatively insensitive to the range ([M]) for m^. However, Fig. [TUl shows that this is not the case for the LET-sensitive 



NGSR N^f^'^ which exhibits strong dependence on mg. Thus we focus our analysis on the fc = NGSR Nlf^' as it 
is less affected by quark-mass uncertainties. 

Our detailed analysis of the non-diagonal GSR begins with an exploration of its consistency with the results of 



(99) 
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Figure 8: Comparison of the complete non-diagonal NGSR N^' (s, r, sq) with that constructed using only the 
leading order chiral contributions (quark condensate and instanton terms) given in (|96p . Central values of the QCD 
input parameters have been employed along with sq = 2.5 GeV^. The upper set of curves are for r — 2 GeV^ and 
the bottom set of curves is for r = 4GeV"'. The two curves overlap almost completely. 
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Figure 9: Comparison of the the NGSR N^^'^' {s,t,so) for the upper and lower ranges of the quark mass specified 
in ([92]) . Central values of the other QCD input parameters have been employed along with t = SGeV'* and 
So — 2.5 GeV^. The two curves overlap almost completely. 



the diagonal cases. In the double narrow resonance model, the non-diagonal NGSR has the form 
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'Attt 
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(s - mlf 
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flgflq + f2gf2q flgflq + f2gf2q 

From the analysis of the diagonal cases, we have found mi w 1 GeV, TO2 ~ 1-4 GeV, and 

/i/ Jgg) ^ 1 _ 0^303 ^ 



Sgg) 
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flg^ + 



flq 



flq + f2q 
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2 ' '2 



1 - 0.593 



/iff^ + f2g^ 
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(97) 
(98) 

(99) 

(100) 



Thus the parameterization of the mixed gluonic-gg system has four couplings of the states to the various currents as 
in Ref. [5J. The four equations ([99]) and (|100p representing the diagonal results determine the four couplings up to 
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Fiffure 10: Comparison of the the NGSR N^J^' (s, t, sq) for the upper and lower ranges of the quark mass specified 
in (jM]) . Central values of the other QCD input parameters have been employed along with r = SGeV* and 
So = 2.5 GeV^. The solid and dotted curves respectively correspond to the lower and upper bound on the quark 
mass. 



an overall sign, leading to two possible solutions for the non-diagonal case 

. (101) 

^ [-3.899 

Apart from the ambiguity arising from the sign of the couplings, all the phenomenological parameters in the non- 
diagonal NGSR (I97|) are determined except for the continuum sq which can be determined by performing a least- 
squares fit of the J, T dependence of (|97p in the region range — 4GeV^ < s < 8GeV^ and 2GeV^ < r < 4GeV'*. 
The best fit for the two cases in (jlOip are shown in Figs. [TT] and [T^ From these Figures we see that the positive case 
in poip is demonstrably most consistent with the QCD prediction. 
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Figure 11: Comparison of the best fit of theoretical prediction for the normalized GSR Nq (s, t, sq) to the double 
narrow resonance phenomenological model. Resonance parameters resulting from the analyses of the diagonal NGSRs 
have been employed (the solution ri = 0.443 from (jlOip . along with mi = 1 GeV and m2 = 1.4 GeV). The optimized 
value of the continuum for these (inputted) resonance parameters is sq = 2.8 GeV^. The upper set of curves are for 
T = 2 GeV* and the bottom set of curves is for r = 4 GeV*. The phenomenological and QCD predictions overlap to 
a large extent. 

At this point we can reach an important conclusion: the non-diagonal GSR is consistent with the results of the 
diagonal gluonic GSR analyses, providing strong evidence for a consistent scenario of mixed gluonic-gg states with 
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Figure 12: Comparison of the best fit of theoretical prediction for the normahzed GSR Nq (s, t, sq) to the double 
narrow resonance phenomcnological model. Resonance parameters resulting from the analyses of the diagonal NGSRs 
have been employed (the solution ri = —3.899 from (|10ip . along with mi — 1 GeV and ni2 = 1.4 GeV). The optimized 
value of the continuum for these (inputted) resonance parameters is sq = 5.2 GeV^. The solid and dashed curves 
represent N^^'^^ (s, t, so) and respectively correspond to r = 2 GeV^ and r = 4 GeV*. The dotted and dashed-dotted 
curves represent the double resonance model and respectively correspond to r = 2 GeV* and r = 4GeV^. 



masses of rrii « 1 GeV and m2 « 1.4 GeV with the heavier state having a slightly larger gluonium content. A 
significant feature of our analysis is the nearly-identical masses that have resulted independently from the diagonal 
gluonic and diagonal qq cases. Although we have demonstrated that the non-diagonal case is consistent with these 
results, an independent analysis of the non-diagonal case would provide further confirmation of the ghionic-qq mixing 
scenario. 

A least-squares fit of the s, t dependence of the non-diagonal NGSR ([97|l in the region range — 4GeV^ < J < 
SGeV^ and 2 GeV^ < r < 4GeV'^ results in sq = 3.00 GeV^, mi = 0.869 GeV, TO2 = 1.39 GeV, r^^''^ = 0.439, 
r^^*-' = 0.561. Fig. [Ol shows that the fit between the QCD prediction and phenomcnological model is excellent, 
with no anomalous deviations indicative of effects from extremely wide resonances. The effect of varying the quark 
mass within the range (|92[) is minimal, so the dominant uncertainty in these predictions arise from the instanton 
parameters. One can estimate the effect of the instanton parameters by recognizing that the instanton contribution 
in (|S7|) is proportional to 



aw wy'w exp 







[ Ji (w^) Y2 {w^) + J2 {in^) Vi {w^)] = / (p^l, pV, p'so) . (102) 



Thus if the instanton size is re-scaled by p — > p/^, its contribution to the NGSR is invariant under the scalings 
So — >■ ^^so, s — > ^^s, r — > ^'*T. Now if these scaling laws are also applied to the double narrow resonance model on 
the right-hand side of Il97p . then it will also be invariant if the resonance parameters scale as mi — > ^^mi, TO2 — *■ ^^W2, 



^(99) _^ ^im) ^ ^im) _j, [To| xhus the results of the best fit analysis will be approximately insensitive to changes in 
the instanton density ric, the resonance strength will be relatively stable under variations in the instanton size, and 
the masses will scale proportionally with the instanton size. This implies that the 15% uncertainty in the instanton 
size will result in a systematic ~ 15% uncertainty {i.e., each mass is shifted in the same direction, leading to the 
same 15% uncertainty in the mass splitting toi — TO2) in the mass predictions. It thus seems reasonable to conclude, 
taking other sources of theoretical uncertainty into account, that the pattern of uncertainty for the non-diagonal case 
is similar to that of the diagonal gluonic case; a 0.2 GeV correlated uncertainty in the mass predictions that leads to 
a relatively stable mass splitting that can deviate by 10% from the central value toi — m2 ~ 0.5 GeV. 

The results obtained from the diagonal and non-diagonal cases are summarized in Table [TJ Taking into account 
the uncertainty in the fitted values of the resonance parameters for the diagonal and non-diagonal NGSRs, we see 
that the non-diagonal case leads to predictions that are consistent with those of the diagonal analyses, so that 



^"Note that overall factors, such as the normalization of the double resonance model and the instanton density, do not modify the 
NGSR and hence do not alter the scaling argument. 
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Figure 13: Comparison of the theoretical prediction for the non- diagonal NGSR N^^"^^ (s, r, sq) with the double 
narrow resonance phenomenological model for the optimized value of the continuum sq and resonance parameters 
(mi, 7712, ri). The upper set of curves are for t — 2GeV'' and the bottom set of curves is for r — 4GeV''. 

all possible GSRs of gluonic and qq currents independently confirm the existence of states of mass ~ 1 GeV and 
~ 1.4 GeV which are mixtures of gluonium and qq mesons. In particular, there is excellent agreement between the 
value of r[^'^^ = 0.439 obtained from analysis of the non-diagonal NGSR and the positive solution r[^'^^ = 0.443 in 
Eq. (jlOip resulting from the diagonal NGSRs. 



Sum- Rule 


TOi (GcV) 


TO2 (GeV) 


ri 


r2 


so (GcV^) 


diagonal: gluonic-giuonic 


0.951 


1.41 


0.303 


0.697 


2.30 


diagonal: qq-qq 


0.943 


1.41 


0.593 


0.407 


2.35 


non-diagonal: gluonic-gf; 


0.869 


1.39 


0.439 


0.561 


3.00 



Table 1 : Analysis results from the diagonal and non-diagonal NGSRs of gluonic and qq currents in the double narrow 
resonance model. Central values of the QCD input parameters have been employed. 

We now examine the pattern of mixing of these states as contained in the quantities rt for the various cases. First, 
we note that although the fitted value of rj^*'' is approximately insensitive to variations in the instanton size, the 
dominance of non-perturbative effects in the non-diagonal NGSR implies that there will be systematic theoretical 
uncertainties associated with truncation of the operator-product expansion and the single-instanton and instanton 
liquid approximation. An uncertainty in rj^*'' of approximately 10%, at the same level as that found in the diagonal 
gluonic case, is therefore appropriate. Thus the value rj^'''' = 0.439 found from the non-diagonal case is in excellent 
agreement with the Eq (|10ip solution r^^*'' = 0.443 that emerges from the diagonal analyses. 

If the couplings obey a single-angle mixing pattern, then one would expect r^^^'' = cos^ 9 — r['''^\ Although such 
a scenario could be possible given the 10% uncertainty in these quantities in Table [l] the non-diagonal case provides 
a more interesting test because single-angle mixing leads to 

(103) 

In this situation, the integral of the right-hand side of (|103p is zero which then requires Mq^'-* (t, sq) = in (|^. 
A value of sq = 4.35 GeV^ can be found to satisfy this constraint over the considered range 2 GeV* < r < 4 GeV*. 
However, as shown in Fig. [T?] the s, r dependence of the QCD prediction Gq^'^'' (s,r, sq) is not consistent with 
mass scales mi w 1 GeV and m2 ~ 1.4 GeV. We thus conclude that the pattern of mixing for the couplings is 
not well-described by a single mixing angle, and hence the situation must be similar to the two-angle scenario that 
has been found for the 77-77' system in the singiet-octet basis [45 . Implicitly this is the same result found in jS], 
where four independent couplings are found necessary in the study of the mixed gluonic-qq system rather than the 
three-parameter system of two couplings and one mixing angle. 




4''^ (S,T,50) 



sin 29 



exp 



(s — mf)'' 



exp 



,2\2 



4t 



20 




Figure 14: Comparison of the best fit of theoretical prediction for the GSR Gq^''' (s, r, sq) to the double narrow 
resonance phenomenological model for a single mixing-angle. Resonance masses mi = 1 GeV and ni2 — 1.4 GeV 
resulting from the analyses of the diagonal NGSRs have been employed. The optimized value of the continuum 
leading to Mq = for these resonance parameters is sq = 4.35 GeV^. The solid and dotted curves represent 
Gq^*^ (s,t, So) and respectively correspond to r = 2 GcV^ and t = 4GeV^. The dashed and dashed-dotted curves 
represent the single mixing- angle double resonance model and respectively correspond to r = 2 GeV^ and r = 4 GeV''. 

In the two-angle scenario the state couplings satisfy [45] 

f[''^ = fg cos Og , /i^^) - /, Sin Og , (104) 
=-fg Sin Og , = /, COS Og , (105) 

which imply that 

ri""^ ^ sinHg , ri'"''^ ^ cosHg (107) 
(gq) _ sin i^9 ~ ^q) - sin jOg + Og) (gq) _ sin (6>g - 6g) + sin {Og -f Og) 

1 2sin(0g-0,) ' 2 2sin(0g-0q) ' ^ ' 

If we use the Table [T] values for the diagonal cases to determine the magnitudes of 9g and 9g and the non-diagonal 
case to determine the relative signs of the angles, then we find the central values 9g = 57° and Og = — 50°li3 This 
solution suggests that there is effectively a single scale for the mixing angle corresponding to a small angle solution 
for 9g + 9g. Following Ref . [T^ , we then define an effective mixing angle a 

Taking into account the uncertainties in the Table [1] values and requiring that solutions to Eqs. p06)) - p08ll be 
consistent leads to Og + 9q — 6.5° ± 1° and a — 54° ± 4°. Thus we find that the effective mixing angle corresponds 
to nearly maximal mixing {a — 45°) where each of the two states are equally coupled to the qq and gluonic currents. 
The deviation of our effective mixing angle from the maximal angle indicates that the heavier (1.4 GeV) state |2) is 
somewhat more gluonic in comparison to the lighter (1.0 GeV) state |1). 



5 Discussion and Conclusions 

Gaussian QCD sum-rules are able to probe hadronic spectral functions over a broad range of energy, and are thus 
ideally suited to exploring the possibility that mixed giuonium-gg states exist amongst the light scalar mesons. We 
have studied the NGSRs for all possible combinations of scalar gluonic and scalar / = (non-strange) qq currents 

There remains an overall sign ambiguity in the mixing angles. 
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(diagonal gluonic, diagonal qq, and non-diagonal gg-gluonic) and find that all three cases independently predict the 
existence of a ~ 1 GeV state and a ~ 1.4 GeV state. This is precisely what one would expect from hadronic states 
that are mixtures of gluonium and quark mesons. Given the uncertainties in our mass predictions combined with the 
result that these predictions represent an upper bound for very broad resonances, it is not clear whether our lighter 
state should be interpreted as the /o(980) or a (at the heavier end of its range pQ) and it is not clear whether the 
heavier state should be interpreted as the /o(1370) or the /o(1500). However, because the approximate 0.5 GeV mass 
splitting between the states is relatively stable under QCD uncertainties, our results do suggest identifying either 
the fighter pair of states [a and /o(1370)] or the heavier pair [/o(980) and /o(1500)] as mixed gluonium-gg states. 

The non-diagonal sum-rule provides important insights into the mixing of the qq and gluonic aspects of these 
two states. Because of chiral suppression factors associated with the light (non-strange) quarks, perturbative effects 
are unable to generate any significant amount of mixing. However, the contributions of the quark condensate 
and instanton effects do not suffer from this chiral suppression and provide the dominant contribution to the non- 
diagonal correlation function, implying that mixing of gluonic and qq degrees of freedom has a non-perturbative 
origin. Qualitatively, this conclusion is similar to one obtained in [44] in which it was demonstrated that instantons 
can lead to a significant mixing between glueballs and (heavy quark) mesons in the pseudoscalar channel. 

The state couplings that result from the analysis of the various GSRs also provide an additional means to examine 
the self-consistency of the scenario of « 1 GeV and « 1.4 GeV mixed states. In particular, the relative couplings 
between the states in the non-diagonal case is constrained by the relative couplings in the diagonal case. The 
independent prediction of these couplings from the non-diagonal NGSR is found to satisfy this constraint extremely 
well, providing strong evidence for the validity of the mixing scenario. 

The state couplings also provide a means to study the pattern of mixing associated with the gluonic and qq 
currents. The resulting pattern is similar to the two-angle mixing that occurs in the rj-rj' system in the singlet-octet 
basis I45j , and results in an effective mixing angle of a « 54°. Because this mixing angle is in the region near maximal 
mixing (a — 45°), there is only a slight preference for the heavier 1.4 GeV state to couple to gluonic channels and a 
concomitantly slight preference for the lighter 1.0 GeV state to couple to qq channels. Indeed, the existence of such 
strong mixing implies that qualitative features that would distinguish pure gluonic and qq states would be obscured 
for strongly-mixed states and the experimental signal of gluonium would thus be subtle. 

In summary, our results provide strong evidence to support the scenario where the mixing of qq and gluonium is 
manifested in the scalar hadronic spectrum as a lighter state on the order of 1 GeV and a heavier state on the order 
of 1.5 GeV m m ini m m Uni HH- in particular, our conclusion that there exists a strong mixing between gluonium 
and qq states is similar to the results of Refs. [ini [H] and the heavier state's slight preference for gluonic channels 
has implications for the gluonium content of the /o(1500) [JS]. 
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